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Abstract 

Material failure by crack propagation essentially involves a concentration of large displacement-gradients 
near a crack's tip, even at scales where no irreversible deformation and energy dissipation occurs. This 
physical situation provides the motivation for a systematic gradient expansion of general nonlinear elastic 
constitutive laws that goes beyond the first order displacement-gradient expansion that is the basis for linear 
elastic fracture mechanics (LEFM). A weakly nonlinear fracture mechanics theory was recently developed by 
considering displacement-gradients up to second order. The theory predicts that, at scales within a dynamic 
lengthscale £ from a crack's tip, significant logr displacements and 1/r displacement-gradient contributions 
arise. Whereas in LEFM the 1/r singularity generates an unbalanced force and must be discarded, we show 
that this singularity not only exists but is necessary in the weakly nonlinear theory. The theory generates 
no spurious forces and is consistent with the notion of the autonomy of the near-tip nonlinear region. The 
J-integral in the weakly nonlinear theory is also shown to be path- independent, taking the same value as the 
linear elastic J-integral. Thus, the weakly nonlinear theory retains the key tenets of fracture mechanics, while 
providing excellent quantitative agreement with measurements near the tip of single propagating cracks. As 
£ is consistent with lengthscales that appear in crack tip instabilities, we suggest that this theory may serve 
as a promising starting point for resolving open questions in fracture dynamics. 
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1. Introduction 



The dynamics of propagating cracks, and their accompanying instabilities ( Fineberg and Marder . 19991 : 
Livne et al.l . l2007l ). remain a fundamental and challenging scientific problem. The most developed theoret ical 
approach to fracture dynamics. Linear Elastic Fracture Mechanics (LEFM) (jFreund . 1998 : Broberg, 199i^ . is 
based on the assumption that deviations from a linear elastic constitutive material behavior takes place only 
in a small region near the tip of a crack. The major prediction of this theoretical framework is that outside 
of this small nonlinear near-tip region, but not very far from it, the displacement-gradient (and stress) fields 
are characterized by a l/i/r singularity, where r is measured from the crack's tip. The intensity of the 
singularity is quantified by the stress intensity factor, that appears as a pre-factor in the singular fields 
and depends on the applied loadings and geometric configuration in a given problem. A corollary of this 
basic result is that the crack tip opening displacement (CTOD) is predicted to be parabolic, where the 
curvature of the parabola is a function of K. 



A central concep t in LEFM is that of the autonomy of the near-tip nonlinear zone (jRicd . Il968at iFreundl . 

1998 HBrober"5ll999l) . The basic idea is that the mechanical state within the near-tip nonlinear zone, which 
is surrounded by the l/^/r singular fields (the "K-fields"), is determined uniquely by the value of K, but 
is otherwise independent of the applied loadings and the geometric configuration in a given problem. This 
implies, for example, that systems with the same K, but with different applied loadings and geometric 
configurations, will be in the same mechanical state within the near-tip nonlinear zone. 

One of the most useful applications of the 1 / ^/F singular fields is the calculation of the energy flow into 
the nonlinear near-tip zone during crack propagation. It is calculated using a path-independent integral, the 
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J-integral (Rice, 1968bl : Freund . 1998 ). that is evaluated within the range of dominance of the l/-y/r singular 
fields. This energy flux, which is quadratic in K , is eventually dissipated at smaller scales near the tip of 
the crack. 

The '^1 \/r singular fields, however, imply that at some distance from the crack's tip the small displacement- 
gradient assumption of LEFM must break down when higher order displacement-gradients become non- 
negligible. This breakdown of LEFM is expected to take place prior to the intervention of irreversible 
deformation and energy dissipation, and thus should be studied in the framework of nonlinear elasticity. 
This physical situation motivates a controlled expansion of a general nonlinear elastic material constitutive 
law in powers of displacement-gradients, beyond the common small displacement- gradient approximation 
of LEFM. Such an approach was taken recently, when a weakly nonlinear fracture mechanics theory was 
developed by deriving the leading no nlinear elastic correction, up to second orde r displacement-gradients, 
to the 1/ y/r singular fields of LEFM I Livne et al.l . l2008t iBouchbinder et al. . 20081) . It is important to note 
that this approach is entirely general as any material exhibits nonlinear elastic response to some degree. 
Moreover, it entails no new constitutive assumptions and/or the introduction of a new intrinsic material 
lengthscale. 

The new theory was te sted against direct measurements of the deformation near the tip of propagating 
cracks ( Livne et al. . 2008h and was sh own to be in excellent quantitative agreement with these experi- 



mental data (jBouchbinder et al.l . 120081 ) . The weakly nonlinear correction to LEFM includes displacement 



contributions proportional to logr and displacement-gradient contributions proportional to 1/r. The logr 
contributions most clearly manifest themselves by a strong departure from a parabolic CTOD in the near 
vicinity of the crack tip. The 1/r strain contributions immediately imply that within the spatial range of 
validity of the weakly nonlinear solution, the displacement-gradient spatial variation is "more divergent" 
than the common singularity. 

The 1/r displacement-gradient singularity must be discarded i n the framework of LEFM, since it gen- 
erates an unbalanced force acting on a line encircling a crack's tip (|Ricd . l1974). In this work we show that 
in the framework of the weakly nonlinear theory this singularity generates no such spurious force and is 
both necessary and physically acceptable. Furthermore, the condition that such a spurious force does not 
exist can be used to retain a key feature of LEFM, namely that of the autonomy of the near-tip nonlinear 
zone, implying that the weakly nonlinear solution depends only on the stress intensity factor K . Finally, we 
calculate the weakly nonlinear J-integral, i.e. we expand the general nonlinear J-integral up to the leading 
nonlinear correction to LEFM, and show that it takes the same value as its LEFM counterpart. This result 
explicitly demonstrates the path-independence of the weakly nonlinear J-integral and, together with the 
notion of autonomy, shows that the weakly nonlinear theory is consistent with all of the common tenets of 
fracture mechanics. 

The weakly nonlinear solution implies the existence of a lengthscale £ that is associated with the break- 
down of LEFM near the tip of a propagating crack, i is defined as the lengthscale at which the 1 /r singular 
terms become non-negligible compared to the dor ninan t 1/ \/ r linear cla stic singul ar terms and marks the 
onset of deformation - depen dent material behavior ( Gaol . Il996t iBuehler et al.l . l2003t 



Buehler and Gaol . l200d 



Bouchbinder and Lol . 120081 ). We suggest that £ may be intimately related to crack tip instabilities that 



currently remain unexplained within the framework of LEFM. Therefore, the set of encouraging results 
obtained in this work leads us to suggest that the weakly nonlinear theory may be a promising starting 
point for unlocking a plethora of open questions in fracture dynamics. This line of investigation is close in 
spirit to other rece ntly proposed approaches that highlight the role of nonlinear elasticity near the tip of 
propagating cracks (jBuehler et al.l . 120031 : iBuehler and "Gaol . l2006l) . 

The s tructure of this pap er is as follows. In Sect. Hjwe describe the recently derived weakly nonlinear 
solution ( Bouchbinder et al.l . 2008) and include a detailed derivation of the theory. In Sect. [3] we discuss 
the difference between LEFM and the weakly nonlinear theory in relation to the 1/r singularity. We also 
demonstrate the autonomy of the near-tip nonlinear zone. This section culminates in a quasi-static weakly 
nonlinear solution. In Sect. [Hwc compare the resulting solution with direct measurements of the deformation 
near the tip of a crack propagating in a neo-Hookean material and demonstrate excellent agreement with 
experiments. In Sect. [Hwe show that the weakly nonlinear J-integral is path- independent, its value coinciding 
with the well-known LEFM result. Section [6] offers some concluding remarks. 
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2. The details of the weakly nonlinear solution 



As noted in the introductio n, nonlinear materi al response at the large strains near a crack's tip in general 
and the experimental results of iLivne et al. I (|2008h in particular, motivate us to formulate a nonlinear elastic 
dynamic fracture problem under plane stress conditions. In this section we derive, in det ail, the weakly 
nonlinear theory, which was previously outlined in the results of iBouchbinder et alj (|2008r ). Consider the 
deformation field (p, which is assumed to be a continuous, difFerentiable and invertible mapping between a 
reference configuration x and a deformed configuration x' such that 

x' = 4>{x) = X + u{x) , (1) 

where u{x) is the displacement field. The deformation gradient tensor F is defined as 

F - V0 , (2) 

or explicitly as Fij = Sij + djUi, where i,j denote Cartesian components. The first Piola-Kirchhoff stress 
tensor s, that is work- conjugate to the deformation gradient F, is given as 

where U[F) is th e strain energy in the deformed configuration per unit volume in the reference configuration 
(iHolzapfeifioOoh . The momentum balance equation is 

V • s = pdttcf> , (4) 

where p is the reference mass density. Under steady-state propagation conditions we expect all of the fields 
to depend on x and t through the combination x — vt and therefore dt — ~vdx- Here x is the propagation 
direction and y is the loading direction. The polar coordinate system (r, 9) that moves with the crack tip is 
related to the rest frame by r= {x — vt)"^ -I- and 9 = tSia~^[y / {x — vt)]. Thus, the traction-free boundary 
conditions on the crack's faces are 

Sxy{r,9 = ±Ti)^Syy{r,9^±TT) . (5) 

We now perform a controlled expansion of the form 

M(r, 9) ~ (r, 61) + e^tt^^) (r, 61) + ©(e^) , (6) 

where e is a measure of the magnitude of displacement-gradients. The expansion is perturbative in the 
sense that 0{t^) contributions to the deformation fields are neglected and asymptotic in the sense that we 
consider a region near the crack tip where C(e^) contributions are non-negligible compared to the leading 
0{e) contributions in this region. Equation ^ implies that we consider first (linear) and second (quadratic) 
orders of elasticity. The order e problem is the standard LEFM one 

/iV^M^i) + (A + /i)V(V • it(i') = pii^i) , (7) 

with the traction- free boundary conditions on the crack's faces at 6* = ±7r 

(A + 2/i)^^-f A5,uW =0 . (8) 

Here A is the plane-stress first Lame coefficient and fi is the shear m odulus. The n ear-tip asymptotic solution 
under Mode I steady-state propagation conditions is known to be ( Freundl . ll998[ ) 



4^V27r 4^(A + fi) 

eu^^\r,e;v) = ' )_ 17^ 6>;« = . (9) 

y ' ' ApL^ ' 4/i(A + Ai) ^ ^ 
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Here Ki is the Mo de I stre s s inte nsity factor (denoted in the introduction by K) and T is a constant known 

1998f) . Note that these parameters cannot be determined by t he asyrn p 

411 



as the "T-stress" (jFreundl . Il998f) . Note that these parameters cannot be determined by t he asyrnptotic 
analysis as they depend on the global crack problem. r2(0;w) is a known universal function ( Freundl . 119981: 



Bouchbinder et al.l . 120081 ) . e in Eq 



can be now defined explicitly as 



4:tiy/2n£{v) 



(10) 



where £{v) is a velocity-dependent length-scale. £{v) defines the scale where order contributions become 
non-negligible compared to order e contributions, which is a property of the nonlinear elastic constitutive 
behavior of a given material. This is a "dynamic" (as opposed to "intrinsic") length-scale that marks the 
onset of deviations from a linear elastic constitutive behavior. 



The order Mode I problem has the following form 



/xV2M(2) + (A + /i)V(V-M 



pu 



(2) 



(11) 



To this order the traction- free boundary conditions at 9 — ±t: now become 

a (2) 



0, 



Xdr 



,(2) 



-I- 



p£k{v) 



= 0. 



(12) 



Here g(9;v) and k{v) are universal functions that depend on the first and second order elastic moduli. 
Generic quadratic nonlinearities of the form d{du^^'^du'-^^^) result in an effective body-force oc in Eq. 
(fTTjl . Note that corrections proportional to KiTr~^^^ / fj, were neglected in Eq. (fTTj) . These corrections lead 
to displacement contributions that vary as r^/^ and thus result in a few percent v ariation of the app arent 
Kj with the angle 9^ This may indeed be observed, cf. Figs. [1] and [H see also iLivne et al. (|2008l) and 



Bouchbinder et al.l (|2008l) 



The solution of the order Mode I problem posed by Eqs. pT|) - (fT^ is obtained in two steps. First we 
obtain a particular solution £Y{9;v) of Eq. (fTTj) . This solution is r-independent and results from solving 
a vectorial linear ordinary differential equation in 9 that depends on g{9] v) and thus on the second order 
elastic moduli. For t h e par ticular case of a neo-Hookean material, T{9; v) (for a general v) was obtained 



Bouchbinder et al. I liooi). Since £T{9; v) does not satisfy the boundary conditions of Eq. we have 



to add to it a solution of the homogeneous equation, obtained by omitting the inhomogeneous term in Eq. 
pip . The resulting equation is just the Lame equation of ([7]). The homogeneous solution is selected such 
that the sum of the particular and homogeneous solutions satisfies the boundary conditions of Eq. p2p . 
Denoting the homogeneous solution by £u, we obtain 



ir,9;v) 



(r,( 



Ki 



Ki 



[u^{r,9;v) + T,{9;v)]. 



[uy{r,9;v) + Tyi9;v)]. 



To obtain u we follow a standard procedure ( Freund . 1998; iBroberel . Il999l) . Using the Helmholtz de 
composition 



u = \/lp + V X 



we obtain 



aldli; + = 



(14) 
(15) 

J,, where Cd.s are 



where = tpz [z is a unit vector perpendicular to the xy-plane). Here ^ = 1 — u /c^ 
the dilatationa l and shear wave speeds respectively and (/3, ip are the commonly used displacement potentials 
( Freund . 1998f ). These equations can be rewritten as two Laplace's equations in the complex variables 



X + ladV = rae 



and 
4 



(16) 



which are coupled only on the boundaries. Note that tan^cj ^ ^a^^stan^ and Vd^s — f^^ — (wsin0/cd,s)^. 
The most general solution of Eqs. (fT5| is given as 



^[t, 9) ^ 3?{i?(Cd)}, ^(r, e) = 3{Q(C)} 



(17) 



where R and Q are two analytic functions. The crucial point for our discussion is that the second boundary 
condition of Eq. (fT^ contains a term that varies as r^^ . Therefore, the solution must be one that has the 
property that Vm cx r~^. The solution with this property is given by 



R{U) = Hd log Cd, Q(C.) = log C. 



Using the relations (cf. Eq. (fTi]) ) 



(18) 



(19) 



we obtain 



Uy{r, 9, v) 



A log r 4 



A. 
2 log 

-50, 



Sa, log r 



log 



sin^ ( 



(20) 



When this sol ution is substituted in the boundary conditions of Eq. (|12p . a single relation between A and 
B is obtained ( Bouchbinder et al.l . 2008 ). The remaining free parameter, say B, was previously determined 



Bouchbinder et al.l (|2008l ) by fitting the solution to the experimental data. In the next section we will 



determine this free parameter theoretically. This result shows that no free parameters exist in the asymptotic 
solution in addition to Kj (and possibly T). This is in accord with the concept of the autonomy of the 
nonlinear near-tip region. 



3. The autonomy of the near-tip nonlinear zone 

In this section we directly address the difference between LEFM and the weakly nonlinear theory in 
relation to the existence of the 1/r singularity. We will show how this difference is used to retain the 
autonomy of the near-tip nonlinear zone. Let us now, for concret eness, consider the neo-Hookean plane- 
stress elastic strain energy functional ( Knowles a nd Sternberg, 'l983l) that describes the constitutive behavior 
of the elastomer gels used in the experiments of Livnc et al. (2008) 

U{F) = ^[F,,F,,+\^ -i] , (21) 
where the out of plane stretch A (not to be confused with the plane-stress first Lame coefficient A) is given 

by 

A = det(F)-i . (22) 

The last relation corresponds to the incompressibility condition, which implies a single elastic modulus /i. 
Incompressibility entirely determines the second linear elastic modulus to be A = 2^, which is equivalent to 
the incompressible linear elastic Poisson ration v = 1/2. For later use we write A(r, 9) of Eq. (1^^ explicitly 
as 

A(r, 9) = ^ , „ , \^ ^ ■ (23) 



The equations of motion corresponding to this energy functional, derived using Eqs. ©-(gl), become 

AiV^^o; + -[dgX^dr(l)y~drX^d0(l)y\=pV^dxxlpx, 



r 



mV'0„ + ^[drX'de(kx-deX'dr(^x\=pv''dxx(l)y, (24) 
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where the assumption of steady-state propagation condition was used. The traction-free boundary conditions 
on the crack's faces, cf. Eq. (O, are 



= 0, 



= 



(25) 



In order to derive the functions g{9; v) and k{v) in Eqs. (fTT|l - (fT^ we expand A'^(r, 6), where A is given in 
Eq. (|23p . to second order in e (see Appendix). Substituting this expansion (see Eq. (|44p in the Appendix) 
into Eqs. (|24p - (P5|) and using the first order solution of Eq. ([7]) we obtain explicit expressions for the 
functions g{0;v) and k{v) in Eqs. (fTT|) - p^ . 

To elucidate the mathematical procedure to follow and in order to derive analytical results, we now 
restrict ourselves to the quasi-static limit v 0. In this limit we obtain 



9xiO;v = 0) 



12 
1 



55 cos d + 208 cos(26i) - 27 cos(36')] , 



gy{e;v = 0) ^ — [-79sin6l-|- 208sin(26')-27sin(36l)], 



Using £/(e',0) of Eq 
appearing in Eq. ([13 



<v = 0) - -— . (26) 
, we obtain the components of the particular (inhomogeneous) part of the solution 



= -i|^cos0 + ^cos(2(?)-4'^os(3^)' 
48 15 16 



61 



26 



y{e;0) = --sin0+-sm{29)--sm{30). 



(27) 



The quasi-static limit of the homogeneous part of solution appearing in Eq. is a bit subtle. Both 
the parameters A and B in Eq. ()18p diverge in the limit i; — > according to 



A 
B 



a + 



7c; 



(28) 



A similar divergence is encountered in the asymptotic LEFM result through the appearance of the Rayleigh 
function (iFreimdI . flOoi iBroberd . [Toooh . Therefore, simply substituting w = in Eq. p3p gives a wrong result 
in which there is only a single parameter [A + B) in the solution. To circumvent this problem, we expand 
the homogeneous part of the solution appearing in Eq. (|20p to second order in v /cs 



u^{r,e-0)~{A + B)\ogr 
Uy{r,e]Q) ~ -{A + B)e + 



1 1 \ \ 1? 
-Blogr- -Aain^e- -Bsin^e] —+0(v^), 

2 8 2 J d 



-AO + —A sin(26') + -B sin(20) 
8 16 ^ ^ 4 ^ ' 



(29) 



where we used Cd = 2cs in Eq. (j20p , the appropriate relation for an incompressible material under plane 
stress conditions. Using Eq. ([25]) . the limit ?; ^ now yields 



Uy{r, 9; 0) 



a log r + b [ log + ^ sin^ ( 



(30) 



where 



a + P and 



(31) 



We now focus on determining the parameters a and b in Eq. (I30p . One constraint on these parameters 
is obtained from the boundary conditions of Eq. (|12p . with k{0) of Eq. (|26p. This impUes 



1 

15' 



(32) 



The determination of the second parameter, b, follows from an interesting feature of the 1/r singularity 
of the displacem ent-gradients. As in quasi-static LEFM, the 1/r singular term in a Williams expansion 
( Williamsl . ll957 ) generates a net force transmitted acros s any line su rrounding the crack tip that acts in the 
propagation direction, when b is not properly selected ( Ricd . 1974 ). Since no such force is applied to the 
crack faces/tip and since inertia is negligible in the quasi-static limit, this apparent force is an unbalanced 
one that cannot physically exist. The condition that no net force is transmitted across a line surrounding 
the crack ti p is uniquely satisfied in quasi-static LEFM by eliminating the 1/r strain singularity altogether 
([Ricd . Il974l ). This is not the case in the weakly nonlinear theory. To see this, we calculate the net force / 
transmitted across a line surrounding the crack tip 



(33) 



where the line is chosen to be a circle of radius r whose outward normal is n = (cos 0, sin 0). Note that 
the first Piola-Kirchhoff stress tensor, s, describes forces in the deformed configuration per unit area in the 
reference (un-deformed) configuration. Therefore, the integration in Eq. (|33p over circles in the reference 
configuration, gives the force / on their images in the deformed configuration. 

To evaluate the integrals in Eq. ((33|) . we derive the components of s in Eq. ([3]) for U{F) of Eq. ([2T|) 



fx / [Sj;:^;?!^; -|- Sj^yTly] r (10 ^ 

J —IT 

r 



= fi {dj<pi - X^Cike-jidicpk) , 



(34) 



where is the two-dimensional alternator. More explicitly, in term of the deformation field (p (cf. Eq. 
(dJ), we have 



/i f cos 9dr(j)x — -^'^ sin 9dr(j)y — sin 6 



Ml _ A3 cos 



Syx = II i COS 9 dr(t)y + \ sm9dr(j}x ^ sill' 



r 

Sxy = fJ- [ sin 9dr4>x + A'^ cos 9dr(j)y + cos Q ^^'^^ — \^ sin^ ^^*^^ 



^yy 



n I sin 9dr4>y ^ cos 9dr(f>x + cos Q^^^tl. \^ sinff ^^'^^ 

r r 



(35) 



To expand s in powers of e, we need to expand (p = x -\- u{x) and A'^(r, 9). To this end we first take the 
V ^ limit in the order e solution appearing in Eq. ([9]) and neglect the T-stress term to obtain ( Freundl . 



19981) 



4^)(r,, 



MW(r,( 



U r 



U r 



7 

- cos I - 
3 \ 2 

13 

-r sm I - 
3 \ 2 



i9_ 
2 

W 



(36) 



This, together with the second order solution in Eq. ([50)1 . yields <p to order e^. We then use the expansion 
of X^{r,6) in Eq. to obtain Sy to second order in e. The integrals in Eq. ((55)) are then evaluated 

analytically, yielding 



fx 



oc 



56 
15 



^ + 4a + 36 



To eliminate unbalanced forces, we demand now that f^ in Eq. 
Eq. (1211), yields 

6 = -^. 
45 



(37) 

33]) vanishes. This condition, together with 



(38) 



Therefore, the complete weakly nonlinear quasi-static solution for a neo-Hookean material is entirely deter- 
mined (cf. Eqs. dni), (ED), dSOl) dSH), dSSD and dMD) and reads 



Ux{r, 0) 



Ki 



4^V27r 

2 



7 

- cos I 
.1 V 2 



361 



4/iV27r 



Ki 



4^V27r 

2 



1 52 / 3 

— log(r)- — hog(r) + -sin2(0) 



13 . f0\ . fiO 
_3,n(-j-sm(- 



103 26 3 

cos (6) + — cos (29) cos (30) 

48 ^ 15 ^ ^ 16 ^ ' 



4^V27r 



(39) 



where the only free parameter is Ki, in accord with the concept of the autonomy of the near-tip nonlinear 
zone. The displacement-gradients emerging from this solution still retain a l/r singular contribution. We 
stress that a solution similar to that of Eq. (j39p can be derived for a general strain energy functional, 
retaining all of its generic properties. Here we have focused on a neo-Hookean strain energy functional both 
for concreteness and to enable comparison with the experiments to follow. 

There is a key difference between LEFM and the weakly nonlinear theory in relation to the 1 /r singularity. 
The displacement fields of Eq. ([3D]) are a solution of the (LEFM) Lame equation of ([7]). The boundary 
conditions of Eq. ([5]) are satisfied by setting a~{). However, in LEFM f^ in Eq. (|33p vanishes if and 
only if 6 = 0. This implies that the l/r singularity en ierging from the gradients of the displacement fields 
of Eq. ([5n|) must vanish in the framework of LEFM ( Ricd . Tl973 ) . In contradistinction, the emergence of 
the inhomogeneous terms proportional to g{9;v) and k{v) in Eqs. (fTT]) - (fT2|) . within the framework of the 
weakly nonlinear theory, makes a crucial difference that gives rise to a physically acceptable l/r singular 
contribution. Thus, not only is this singularity physically acceptable, but it is also necessary. 

Finally, we note that Kj in Eq. ([5^ cannot be determined by the asymptotic solution, but rather entails 
either solving the global fracture problem or a direct comparison with an experiment. We will follow the 
latter in the next section. 



4. Comparison to direct measurements 



We now compare the weakly n onlinear solu t ion to the direct measurements of the deformation near a 
propagating crack tip re ported in iLivne et al. ( 2008h . Although a similar comparison was performed in 
Bouchbinder et al. 

' (l2008l) . the resulting agreement of the theory with the experimental data was obtained 
using B of Eq. (fTS]) as a free fitting parameter. We now test whether the solution appearing in Eq. ([5^ , 
where there are no free parameters, still agrees with the experimental observations. 

To perform this comparison , two issues should be addressed. First, the lowest propagation velocity data 
in the experiments reported in iLivne et al.l (|2008r ) correspond to v = 0.2cs, while the solution in Eq. ((39|) 







2 4 
r (mm) 
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Figure 1: (Color online) Inset: The measured Ux{r,0) at v = 0.2cs (circles) fitted to the weakly nonlinear solution (solid line) 
with Ki = 1040Pay'm and T = — 2800Pa. Main panel: The corresponding measurements of eyy{r,0) = dyUy(r,0) (circles) 
compared t o the weakly nonlinear solution (solid line) using the same Kj and T. The LEFM prediction (analysis as in 
iLivne et alJ l|2008l \') is added for comparison (dashed line). 



corresponds to the limit v ^ 0. We therefore replace the quasi-static order e part of Eq. ([39|) . i.e. the 
standard K-field, by its dynamic counterpart for v = 0.2cs (see Eq. Second, in these measurements in 

the same region in space where the order terms are non-negligible compared to order e terms, there exists 
also a non-negligible T-stress contribution ( Bouchbinder et al.l l We, therefore, add the T-stress terms 



appearing in Eq. ([9|) (with A = 2/i to account for incompressibility) to Eq 

The resulting weakly nonlinear solution, with Kj — lQ4QPa,^/rn and T = — 2800Pa, is compared with 
the experimental data for Ux{r,0) and eyy{r,0) — dyUy{r,0) in Fig. [TJ The agreement between the weakly 
nonlinear theory and the experiment is excellent and is a significant improvement to the LEFM result 
(added to Fig. [T] for reference). Moreover, the values of Kj and T used here are close (to within 3% 
and 10%, respectively) to the ones used i n Fig. 1(a) of Bouchbinder et al. ( 2008f l. demonstrating that the 
experimentally selected free parameter in Bouchbinder et al.l ( 20081 ) is consistent with the one determined 
here theoretically. 

For completeness, we compare the CTOD predicted by the weakly nonlinear solution with the measured 
one in Fig. [21 The agreement is favorable and should be compared with Fig. 2(a) in [Bouchbinder et al. 
POOST ). This result explicitly demonstrates the existence of logr displacement terms that are responsible for 
the deviation from a parabolic CTOD observed in Fig. [H We conclude that the weakly nonlinear solution 
in which only Ki (and T in this particular ca se) is a free parame t er agr ees well with the experimental data 



and validates the previously reported results (jBouchbinder et al.l . 1200 



5. The weakly nonlinear J-integral 



We further substantiate the solution in Eq. ([5^ by demonstrating the path-i ndependence of the 
weak l y nonlinear J -integral. The dynamic J-integral for nonlinear elasticity reads (jGurtin and Yatomil . 



1980l : iFreundl . Il998[ ) 



J 



1 



U{F) + -pdtu^dtUi vnx + SijUjdtu 



dC 



(40) 



Here C is a contour that starts at one traction-free face of the crack, surrounds the tip, ends at the opposite 
traction-free face and translates with the crack tip as it moves, n is an outward unit vector on C and U{F) 
i s given i n Eg. ([^ . This J-integral is path-independent for steady-state propagation, i.e. when dt = —vdx 



([Freundl . 119981 ). We note that the J-integral of Eq. is a generalization of the common linear elastic 

J-integral, where U{F) is a general nonlinear strain energy functional, s is the first Piola-Kirchhoff stress 
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-0.25 -0.2 -0.1 -0.05 

(^^(r,7i) (mm) 



Figure 2: (Color online) The measured crack tip profile (0y(r, itTr) vs. tf>x{r,iT)) at v = 0.2cb (circles). Shown are the parabolic 
LEFM best fit (dashed line) and the profile predicted by the weakly nonlinear solution (solid line) with Kj = lOTOPay'm and 
T = — 2800Pa. The values of K[ used here and in Fig. [T]difler by a few percent, which is consistent with the omission of a 
correction of order T/fi mentioned below Eq. H12| l. 



tensor and p is the reference density ( Freundl . 1 19981) . Our aim here is to calculate the integral in Eq. (|40p 
to third order in e, which is termed here the weakly nonlinear J-integral. 

We are interested in the energy release rate G = J/v under quasi-static conditions, i.e. for v ^ 0. 
Therefore we have 



G 



lim — : 



■ [FijFij+X-- 



-3): 



rde 



(41) 



where we choose C to be a circle of radius r ~ £. We then expand and A"^ to order e^, expand Sy in Eq. 
([55)) to order and use Eq. ([M]) to obtain the integrand of Eq. (|4ip to order e^. The integration, which 
involves simple trigonometric functions, can be performed analytically, yielding 



Using e of Eq. Uni) in Eqs. (gUl-dlT]) we obtain 



3fi 



(42) 



(43) 



This is precisely the known LEFM expression (| Freundl . Il998h . This result is expected and shows that the 
weakly nonlinear solution in Eq. (|39p is consistent with the path-independence of the J-integral up to 
order e^. In fact, we must have G*^'^-' — in order for the J-integral to take the same value obtained in the 
asymptotic LEFM region ( Freundl . 1998h . Moreover, the order integrand in Eq. (pijl (which results in G^"^' ) 
is more singular than l/r, which implies an r-dependence that would have violated the path-independence 
property if G'^-* t^O. Thus, again, we must have G^^-* =0, as was obtained explicitly. We conclude that the 
weakly nonlinear solution appearing in Eq. p9p results in a path-independent weakly nonlinear J-integral, 
consistent with the path-independence property of the general J-integral of Eq. (|40p under steady-state 
conditions. 
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6. Concluding remarks 



We demonstrated that the 1/r singularity is a necessary and physically acceptable solution in the frame- 
work of weakly nonlinear fracture mechanics, in contrast with linear elastic fracture mechanics where such 
a singularity does not exist. We showed that the resulting weakly nonlinear solution is compatible with the 
concept of the autonomy of the near-tip nonlinear zone and that the weakly nonlinear J-integral is path- 
independent, taking the same value as its linear elastic counterpart. These results show that the weakly 
nonlinear theory is consistent will all of the tenets of fracture mechanics. 

In the second part of this work we focused on the quasi-static limit, mainly for the sake of illustration 
and in order to derive explicit analytical results. The dynamic weakly nonlinear solution is given explicitly 
in Eqs. (fT^ l) and ([^IB. wh ere the parameters A and B are related by the boundary conditions of Eq. (|12p 
(jBouchbinder et al.l . 2008f l. The remaining parameter, say B, should be determined by a similar procedure 
to the one presented in this work. Note, however, that the appearance of a net force transmitted across a 
line surrounding the crack tip is not necessarily ruled out in the dynamic case {v > 0) since this force can be 
balanced by material inertia (the rate of change of linear momentum) . It is clear though, that the condition 
that the displacement-gradients 1/r si ngularity does not generat e any unbalanced force is a necessary and 
sufficient condition to determine B. In Bouchbinder et air ( 2008h . B was determined as a fittin g parameter 
to the experimental data. Here we showed, for the smallest propagation velocity reported in iLivne et al.l 



( 2008 ): Bouchbinder et al. ( 20081) . that the two proced ures yielded very similar results, giving us confidence 



p roced u 

that the results presented in [Bouchbinder et al. for higher propagation velocities are valid. 



The existence of an implicit lengthscale i{v) associated with the weakly nonlinear theory may open the 
way to understand the origin of crack tip insta bilities, in line with similar ideas that were proposed recently 
( Buehler et al.L 12003 : Buehler and l2006h . There is, in fact, some experimental evidence that points 
in this direction. Livne et al. I (|2007^ reported that, when the micro-branching instability is suppressed in 
gels, sinusoidal oscillations of a crack's path occur beyond a critical velocity of 0.9cs. The wavelength of 
these hig h-velocity crack path osc illations is consistent with the mm-scale £{v) that arises at the relevant 



veloc ity (IBouchbmder et al 



size (Livne et al. 



20081 ). Moreover, this wavelength was sho wn to be independent of the sampl e 

2007f ). as opposed to the LEFM-based prediction (jBouchbinder and Procaccial 120071 ) . 
This observation implies that the wavelength does not originate from the system's geometry, but rather 
from a dynamical effect. In addition, experimental observations of the micro-branching instability in gels 



( Livne et al.l . l200"5h indicated that the typical micro-branch length also seems consistent with £{v) in the 
relevant velocity range. 

This set of observations suggests that the weakly nonlinear theory should be seriously considered as a 
possible route to understand crack tip instabilities. In this regard, we reiterate that the lengthscale £{v) 
marks the breakdown of LEFM and thus the onset of deformation-depende nt material behavior. The latter 
includes, for exarnple, h yperelastic effects such as variable lo cal wave-speeds (iGao , 19961 : Buehler et al. 1. 120031: 
2006f ) and variable local response times ( Bouchbinder and that were suggested 



Buehler and Gao 



as a possible instability mechanisms. The lengthscale £{v) is appealing in this context also from the point of 
view that it is independent of dissipative mechanisms and/or structural lengthscale s . Thi s may be important , 
in view of the existence of identical crack tip instabilities ()Fineberg and Marder . I1999I : iLivne et al.l . l2005^ 
in several amorphous brittle materials with wholly different dissipative mechanisms and micro-structures. 
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A. The expansion of X^{r,6) to second order in e 

In this Appendix we present the expansion of X^{r, 6), where A(r, 6) is given in Eq. 



to second order 



1 - 



3 ( cos 09,41) - ^-^dgum + 3 (sin 09,41) + ^^e^) 



3 { cosOdrU^^ ~ - 3 (sin09,42) ^ S^a,^^) 



+ 6 (cos^ 0(9,41) )2 +sin2 0(9,,4i))2 + £^(a,4i))2 + !!^(a,4i))2 



sm 



(44) 



+ 6 ( sin(20)9,«ii)9,^W + 5^9,i.W9,^W - ^^S^utDS^^ziD - '^^deu^^^ deu\^^ 
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^9,4^)9,41) - ^ 



■9e4')a,4i) ) + 3 ( i9e4')5,MW - -dru'^^^dg^^'^ 



0(e3 
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